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ABSTRACT 

Orbits in different dispersive billiard systems, e.g. the 3 disk system, are mapped 
into a topological well ordered symbol plane and it is showed that forbidden and 
allowed orbits are separated by a monotone pruning front. The pruning front can 
be approximated by a sequence of finite symbolic dynamics grammars. 
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1 Introduction 



The aim of this article is to find an effective description of all allowed orbits in a 
classical billiard system. The systems considered are a point particle on a plane, 
moving freely and reflecting elastically off dispersing walls. In most of our exam- 
ples the walls are the borders to circular shaped reflecting disks. Such systems have 
recently been in focus of much theoretical interest, because of their utility in stud- 
ies of the quantum mechanics of classical chaotic systems H> H> H. Semiclassical 
calculations of the quantal spectrum of classically chaotic systems require summa- 
tions over all periodic orbits of the classical system®, and it is essential to have an 
effective description of the allowed classical orbits. 

The solution we propose here is to apply techniques developed in the study of 
one dimensional unimodal maps, and of two dimensional Henon-type maps. In the 
unimodal maps the ordering of allowed orbits along the parameter axis is given by 

the MSS theory 0, i.e. the kneading sequence of the critical point 0. From the 
kneading sequence one constructs a binary number 7 C and the interval (7 C , 1) is the 
primary pruned region of the unimodal map. An orbit is forbidden if the orbit, 
or any of its images or preimages yields a binary number 7 in the pruned region 
(7c 1). A convenient visualization of the allowed orbits is afforded by a binary tree 
on which every allowed orbit corresponds to a path down the tree. The forbidden 
orbits correspond to the forbidden branches, i.e. the branches that have to be 

pruned from the tree of allowed orbits. 

A unimodal map such as the logistic map has a parameter space window for each 
stable periodic orbit. Within such window the kneading sequence is periodic, one 
can find a finite number of forbidden symbol sub-sequences, and construct a new 

unpruned alphabet^. One can also describe allowed sequences as walks on a finite 

Markov graph constructed from the pruned binary tree0. Chaotic dynamics does 

not in general have a finite Markov partition, but we believe that the statistical 
properties can be approximated from calculations on repellers corresponding to 
converging sequences of nearby windows with stable orbits. 

In ref. Cvitanovic, Gunaratne and Procaccia have extended these methods to 



two dimensional maps like the Henon mapt^ and showed that numerically one can 



obtain a connected primary pruned region in a two dimensional space (7, 5) con- 
structed from the properly ordered symbolic future and the past of the orbits; the 
border of the area is the pruning front. The conjecture is that this curve is mono- 
tone, and that there are no other forbidden orbits than those that pass through 
the primary pruning region. One can approximate the grammar of forbidden orbits 
by forbidding finite sequences that correspond to rectangles in the symbol plane. 
Each such approximation gives a finite number of finite length forbidden symbol 
sequences. From these sequences one can make successively refined Markov par- 
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titions and obtain well converging estimates of averages such as the topological 
HT2l ITU H4i 

entropy" — !' 1 — !' 1 — 11 . The pruning front is the symbol plane representation of the 
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primary homoclinic tangencies introduced by Grassberger and KantzU — U in order 



to partition the Henon map non-wandering set into a "left" and a "right" side. 
1.1 Symbols in billiard systems 

In billiard systems the choice of symbols is much easier; a symbol for bounce 
off each concave wall, e.g. a disk or a branch of the hyperbola, is natural and 
unambiguous E3> ^ 0' 0. Billiards considered here are hyperbolic systems with 
only unstable orbits; the stable and unstable manifolds of the orbits are never tan- 
gential to each other. An orbit in the billiard is unstable also at the point where it 
is pruned, in contrast to the Henon map and smooth Hamiltonian potentials, where 



pruning always involves inverse bifurcations and stable orbits H' E3. 

Construction of billiard symbolic dynamics proceeds in 2 steps: 1. define 
a symbol plane that is topologically well ordered; 2. determine the boundary 
orbits between the legal and the illegal orbits. The well ordered symbol plane is 
defined by introducing a new symbol alphabet and associating with every orbit a 
pair of coordinates constructed from the alphabet. A similar well ordered symbolic 
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alphabet has been introduced for a different scattering problem by Troll The 



second issue is simple. We show that for the billiards one starts a tangential bounce 
from every point on the wall and these are the boundary orbits. 

The symbol plane we construct is complete for a system with a finite number 
of hard disks situated far from each other in the configuration space. A physically 

motivated example is a particle bouncing between 4 hyperbolas^' ^3. This billiard 
do not have a obvious limit of a complete Cantor set but we use the same symbol 
plane as the 4-disks billiard and orbits are pruned the same way. 



2 Covering alphabet for dispersing billiards 

For simplicity we choose to describe the covering alphabets for dispersing billiards 
by using billiards with a finite number of identical circularly shaped walls. In general 
the walls do not have to be circular, but they do have to be dispersing. When they 
are identical, the configuration space symmetry enables us to reduce the number of 
symbols and to simplify the problem. 

If we take a billiard system with sufficiently separated disks in the plane there 
exist an orbit which visits the disks in any desired order. 

Let the disks be enumerated by s G {1,2,..., iV} and denote the bounce of the 
point particle on disk s at the discrete time t by symbol s t . Dispersing (concave) 
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walls can not have two succeeding bounces off the same wall, i.e. the subsequence 
_ss_ is forbidden, but all other sequences are allowed for sufficiently separated disks. 

Given the infinite symbolic past of the particle bouncing on disk so at time 
t — 0, that is . . . s_3S_2'5-i, and the infinite symbolic future . . ., the position 

of the bounce x and the outgoing angle measured from the normal vector of the 
wall are uniquely determined. The phase space coordinates we use are the position 
and the angle (2, 0) of the bounce off disk s . 

If the N-disk system is open, then almost all points (x, 0) correspond to a particle 
that has entered the system from outside at some earlier time and is going to escape 
from the system at some later time. The union of all points contained within the 
billiard for all times is the non- wandering set A. For a sufficiently separated disk 
billiard, A is a Cantor set in the phase space similar to the non-wandering set for 

the Smale horseshoe Each point in this Cantor set is described by a unique 



symbolic past and future. However, if we use the disk symbols s t to describe a 
point in the Cantor set, there is no obvious ordering of the symbols from "small" 
symbols to "large" symbols that would reflect an increase of the position or angle 
coordinate in the phase space. A symbolic description with a natural ordering of 
this kind we call well ordered symbols and we show how these can be constructed 
by a few examples. Other configurations of disks can be worked out in a similar 
way. 



2.1 3 disks symbolic dynamics 

Let the billiard system be 3 equal disks with radius 1 and the distance r between 
the centers, enumerated anticlockwise as in figure [I]. Without loss of generality we 
assume that the particle bounces off disk number 1 at time t = (s = !)• The 
position x is the angle describing the bounce on disc 1 and is the angle from 
the normal vector. Define A^~ as the union of points in phase space that at least 
give t bounces in the billiard after bouncing off disc 1. Thus Af is all (x, 0) of 
disk 1 corresponding to a particle that bounces either off disk 2 or off disk 3 at 
time t = 1. A~£ consists of two diagonal strips in phase space, with decreasing 
when x is increasing. The lower strip in the phase space is the union of all points 
in phase space corresponding to a particle bouncing off disk 2 at time t — 1 and 
the symbol sequence for this strip is SqSi = 12. The upper strip is the union of all 
points corresponding to a particle that bounces off disk 3 at time t — 1, and this 
strip has symbol sequence SqSi = 13. The order of the two strips with increasing x 
or is So^i = {12, 13}. A^ is four diagonal strips, two inside each of the two strips 
of Af. The symbol description of the four strips ordered after increasing values of 
x or is s sis 2 = {121, 123, 132, 131}. Table [l] shows the symbol sequences for 4 
bounces ordered after increasing values of x and 0. The limit of Af when t — > 00 
is a Cantor set of lines with a unique labeling S0S1S2 . . . for each line. 
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Let new symbols be denned from a combination of two symbols s t ~i and s t , 
and such that the ordering of the Cantor set lines in the phase space is preserved. 
In table [I] new symbols w t G {0, 1} are written together with the old symbols for 
the first 4 bounces. ^From the symbols W1W2W3 ... a rational binary number 7 = 
O.W1W2W3 . . . = Y^t^i Wt/2 t is constructed. This symbolic coordinate 7 is increasing 
with increasing values of x and increasing values of <p. To construct the symbols Wt 
let first v t denote if two consecutive bounces s t -i and s t take place in a clockwise or 
anticlockwise direction. The symbols v t are not well ordered because a bounce in a 
dispersing (concave) wall reverses the ordering in the phase space. At each bounce 
one has to invert the symbols v t and this gives the algorithm: 

V t = S t - Sf-l 

if v t < 1 then v t = v t + 3 , . 

f v t -l if todd W 

Wt = \ n -r , 

y 2 — v t 11 t even 

The new symbols are constructed to reflect the ordering of the lines A+ in the 
phase space. Exchanging odd and even in algorithm ([!]) gives symbols w' t = 1 — Wt, 
and a binary number 7' = O.w'^^w^ . . . decreasing with x and 9 and therefore also 
the symbols w' t are well ordered. ^From a orbit described by symbols St we have 
two mappings to the two well ordered symbols Wt and w' t . 

The particle bouncing at (x, 0) also has a past, and a symbol sequence describing 
this past. Let A^T be the union of points in phase space corresponding to a particle 
that has at least t bounces in the billiard before arriving at (x, <p). Then A^f is the 
union of all points arriving at disk 1 after being bounced off disk 2 or off disk 3, and 
Ai is two diagonal strips in the phase space. The incoming angle has opposite sign 
to the outgoing angle and the strips in Aj" have increasing with increasing x. We 
get the new symbols Wt from algorithm (jl]) with t < 0. The symbolic coordinate for 
the past is 5 = • • • = J2t^i Wi-t/2*. The coordinate 5 is increasing with 

increasing value of x but S is decreasing with increasing value of 0. If we choose 
the value 7' for the future we have 5' = O.w'qw'^w'^ ... for the past. 

Let Af = A^ H A^. Then A 2 consist of 16 = (2 • 2) 2 areas in the phase space. 
Each of the 16 disjunct sets of A2 has a unique enumeration by the 4 binary symbols 
{w-i, Wo, Wi, 102}- The set A = Aoo is a Cantor set and each point in the set is 
represented by a bi-infinite symbol string . . . u7_ 2 ^-iu ; ow ; iU'2U'3 .... The set A is the 
non-wandering set of the billiard. 

The symbolic plane is the unit square (7, 6); < 7, 5 < 1. Each non-wandering 
orbit is represented by a point in this plane, and each point in the plane is one of 
the non-wandering orbits. As mentioned earlier, we actually get two points in the 
symbol plane from one orbit. Because of the symmetries of rotation and reflection 
in the billiard, we get a number of different orbits in the billiard from one point 
in the symbol plane, but all the orbits have the same length and stability. The 
symbolic plane is a representation of the phase space of the well separated billiard, 
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with all gaps in the Cantor set removed. This is a very convenient space to work in; 
the curved lines in the phase space become straight lines and this space does not 
change if the distance r changes or if the borders of the disks are slightly changed 
as long as the symmetry and concavity is kept. 



The billiard can be reduced to a fundamental domain I— I. The fundamental 



domain is one sixth of the original 3 disk billiard and it is tiling the whole billiard. In 
our phase space the fundamental domain is the part with x < tt/6 and is constructed 
in the symbol plane as follows: The symbol sequence . . . s_ 2 S-iSoSiS2 . . . gives two 
symbolic coordinates (7, 5) and (7', 5') but if 7 + 5 7^ 1 then one point is in the 
fundamental region and the other point is in the region x > tt/Q. We find the two 
points (7, 5) and (7', 5') and choose (7, 5) if 7 + 5 < 1 or choose (7', 5') if 7 + 5 > 1. If 
7 + 5 = 1 then both (7, 5) and (7', 5') are on the border of the fundamental domain, 
and as our convention we choose (7,5) if 7 > | and choose (7', 5') if 7 < |. If we 
have a billiard without symmetry we only use one map from st to Wt and we do not 
have any fundamental domain. 



2.2 Symbolic dynamics for N disks on a circle 

As a generalization of the 3 disk billiard let N equal disks have the center of each 
disk on a large circle and let the distance between centers of neighbor disks on 
the large circle be r. Then A* is (N — 1) strips in the phase space. The well 
ordered symbols w t G {0, 1,2 ... , (N — 2)} are constructed from the anticlockwise 
enumeration of the disks s t G {1,2,..., N} using the algorithm 

s t - St-i (2) 
v t < 1 then v t = v t + N 
j v t — 1 if t odd 
1 N — v t — 1 if t even 

and the opposite ordered symbols are w' t = N — w t — 2. When N = 3 this is the 

same algorithm as ([I]), and for N = 4 this is the algorithm in ref.ESI. ^From wt we 
construct base (N — 1) symbolic coordinates 7 = O.W1W2W3 . . . = J2t^i Wt/(N — 1)* 
and 5 = O.W0W-1W-.2 ■ ■ ■ = W\-t/{N — 1)* where < 7, 5 < 1, and in a analog 
way 7' and 5'. 

The reduction to the fundamental domain in the symbol plane is the same as 
for 3 disk. If N is even there is one period 2 orbit with only one representation in 
the symbol plane (7,5) = (7', 5') = (1/2, 1/2). 



if 

w t 



2.3 N disks with a center disk 

Let the billiard be a configuration of N disks on a large circle as in the billiard 
above and in addition one disk in the center of this large circle. The radius of the 
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disks are 1 and the distance between neighbor disks are r. The disks on the circle 
is enumerated anticlockwise from 1 to N and a bounce off the disk in the center is 
given the symbol (JV + 1). Figure |2] shows this configuration with JV = 6. If the 
number of disks on the large circle is even, then, because of the disk in the center, 
a point particle can not bounce between two disks opposite to each other on the 
large circle. ^From the (JV + 1) symbols of the disks, we get (JV — 1) well ordered 
symbols, and consists of (JV — 1) strips in the phase space. With N (even) disks 
on the large circle and disk number (JV + 1) in the center of the large circle the 
algorithm defining the well ordered symbols w t G {0, 1,2 ... , (JV — 2)} is 

if s t = (N + l) then w t = (N — 2)/2 

else if s H /(JV + l) then 
v t = s t - S t -l 

if vt < 1 then v t = v t + JV 



1 if todd 

Vf — 1 if t even 



else if s t _i = (JV + 1) 

V t = S t - S t -2 

if v t < -{N -2)/2 then v t = v t + N 
if v t >(N-2)/2 then v t = v t - N 



(3) 



(N-2)/2 + v t if todd 
(JV - 2)/2 - v t if t even 



w t = 

The configuration with N = 6 can be looked at as a first step toward a descrip- 



tion of the Lorentz g triangular lattice with a hard disk in each lattice 



point and a point particle scattering in the lattice. 

If the number of disks iV on the large circle is odd, a point particle can reach 
all other disks after bouncing off one disk when the disks are sufficiently separated. 
The algorithm giving the symbols w t G {0, 1,2 ... , (N — 1)} is 

if St = (N + 1) then w t = (N - l)/2 

else if s t _i ^ (N + 1) then 
v t = s t - S t -! 

if v t < 1 then t> t = v t + N 

if todd 
— 1 if t even 

if w t >{N -l)/2 then iw t = w t + 1 (4) 
else if s t _i = (JV + 1) 

V t = S t - S t -2 

if u t < -(JV - l)/2 then u t = f 4 + JV 
if v t > (N — l)/2 then v t = v t — N 




(N-l)/2 + v t if J odd 
(JV - l)/2 - u t if t even 
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2.4 Orientation exchange at a bounce 



One can explain why all the algorithms reverse the ordering of symbols at each 
bounce by a detailed description of a bounce. Assume two particles move along 
two arbitrarily close parallel lines and bounce off the dispersing (concave) wall at 
an angle 7^ as in fig. |3|. The particle that first hits the wall bounces off the 
wall along a new direction. The second particle crosses the first outgoing line and 
then bounces off the wall. The two outgoing lines cannot cross each other as long 
as the wall is concave or straight. The particle which is on the left side before the 
bounce, moves on the right side after the bounce, and this causes a change in the 
orientation at the bounce. If the angle = 0, the orbits do not cross but because 
the velocity is reversed, the orbit that was on the left side before the bounce is on 
the right side after the bounce. 

For a bounce off a focusing (convex) wall as in fig. |], the two parallel orbits cross 
twice if 7^ 0. There is one crossing after the first particle has bounced off the wall 
and one crossing after the second particle has bounced off the wall. The particle 
which is on the left side before the bounce is on the left side also after the bounce, 
and the orientation does not change. This is true also for = when the orbits 
cross once after the bounces. These bounces take place e.g. for a point particle 
moving in the stadium billiard. Orbits bouncing only off the semicircles par ts o f 

the stadium walls can therefore be described by an ordered rotation numberl23. 
The description of all orbits in the stadium billiard in well ordered symbols are 

given in ref. 



2.5 A condition for sufficiently separated disks 

We want a condition on the geometrical construction of the non-wandering set to 
distinguish between sufficiently separated disks and the case of pruning. We find 
that the disks are sufficiently separated if condition 1 below is true. 

The definition of Af in the general case is as follows: Af consists of M strips 
and each strip m + is the union of points (x, 0) from which a straight line starting 
at x with angle hits a particular disk. As a straight line may go through other 
disks the M strips are not necessarily disjoint. In the same way we define Aj~ as the 
M strips where each strip m~ is the union of points (x, 0) where a line from point 
x with angle —0 hits a particular disk. We call the intersection of a strip m + and 
a strip m~ a rectangle. The set Ai = PI A^f then consists of M 2 rectangles not 
necessarily disjoint. The construction of At with T > 1 follows from demanding 
that the outgoing angle at a bounce is equal to the incoming angle, but allowing 
the straight lines to go through a disc. Then A^ consists of M 2T rectangles. 

Condition 1 There exists a number < T < 00 such that At consists of M 2 
disjoint areas where each area is inside one of the M 2 rectangles of Ai . 



8 



The iteration of the M 2 disjoint areas corresponds to one more bounce and gives 
that each of the M 2 disjoint areas contains M 2 new disjoint areas. The M 4 rectan- 
gles of A 2 then contains M 4 disjoint areas. By induction we find that Kt+t 1 gives 
M 2T disjoint areas inside the M 2T rectangles of Ajv. From this it follows that even 
if the rectangles of Ajv overlap each other, the part of the non-wandering set belong- 
ing to the rectangles do not overlap. A symbol string (w-t>+iW-t>+2 ■ ■ .Wt'-iw' t ), 
with w t G {0, 1, . . . , (M — 1)} uniquely correspond to one rectangle when describing 
the part of the non-wandering set in this rectangle. It then follows that the disks 
are sufficiently separated. 

If there does not exist a T according to the condition it may be that an infinites- 
imal change of the parameter gives a < T < oo and this is the critical parameter 
value where pruning starts. If non of the above is true, then the rectangles always 
overlap and there is pruning in the system and not all symbol strings corresponds 
to an orbit. 

3 Pruning 

When the distance r between the disks is small, a subset of the orbits are forbidden 
(not admissible). We conjecture that only two kinds of forbidden orbits exist in the 
dispersing billiards; orbits passing through a forbidden region (e.g. through a disk) 
and orbits going into a forbidden region and bouncing off the wall from the focusing 
side. Figure |^ a) shows a part of two legal orbits; one passing the dispersing wall 
and one bouncing off the wall on the dispersing side. Figure [5] c) shows a part of 
two forbiddens orbits; one passing through the dispersing wall and the other going 
into the forbidden region and bouncing off the wall from the focusing side. Figure [| 
b) shows that the limit orbit of both orbits is a line tangential to the disk. An orbit 
is therefore pruned together with at least one other orbit. 

The orbits bouncing off the wall from the focusing side are orbits with \<f>\ > 7r/2 
in the phase space. In the phase space for sufficiently separated disks all parts of A 
have angle \<p\ < tt/2. When r decreases, some points in A move closer to the lines 
|0| = 7r/2 and the pruning starts at r c when the points in A with largest value of 
1 4>\ reach these lines. 

When r — > r c with r > r c , the value of T in condition 1 tends to oo. For 
r < r c the different rectangles of A t always have some overlapping. This overlapping 
corresponds to forbidden regions in the symbol plane. The forbidden regions that 
are the simplest to describe are the two regions consisting of points in A with 
|0| > 7r/2. The borders of these regions in phase space are the lines \<f>\ = it/2. In 
symbol plane the two regions correspond to two areas, one in the upper left corner 
and the other in the bottom right corner. The border is the symbolic representation 
of all points in A with <fi = tt/2. To obtain this border, we scan through x values 
with the angle <fi = ±7r/2 and if r > 2 (an open system) in numerical work we only 
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keep points bouncing more than 20 bounces both in the future and in the past. 

As both 7 and 5 increase with the value of x, the border is monotonously 
increasing in the symbol plane. 

The forbidden regions described above contain only the forbidden orbits bounc- 
ing from the focusing side and not the orbits passing through a forbidden region. 
We choose to take the pre-image of this described region as our primary forbidden 
region and the pre-image of the border as our pruning front. The other family of 
forbidden orbits, that is orbits passing through a disk, is an other primary forbidden 
region and has a second pruning front that is very similar to the first pruning front. 
An orbit on this second pruning front has the same symbol sequence as an orbit 
on the first pruning front but without the one symbol for the tangential bounce. 
The two primary forbidden regions have a straight line as a common border and 
this line corresponds to a gap between two bands in Af for sufficiently separated 
disks. We consider the two primary forbidden regions to be bounded by one pruning 
front. This region can be described as the primary overlapping in A. In one of the 
examples we will show that there may be more than one of these primary forbidden 
regions, but all are constructed in the same way. 

The iteration of the points in the symbol plane is a shift operation of the symbols. 
This shift operation is slightly more complicated for well ordered symbols than for 

the original disk symbols, see ref. The primary forbidden region can by the 
shift operation be mapped forward and backward in time. The union of all images 
and preimages of the primary forbidden region is dense and takes the full measure 
in the symbol plane; Their complement, the union of the legal orbits has measure 
zero and is a Cantor set in the symbol plane. 

iFrom the primary forbidden region we can read off all the finite sequences that 
can make an orbit forbidden. If the well ordered alphabet has M symbols, we divide 
the symbol plane into rectangles with side length M~ k . Each rectangle has a unique 
labeling W-k+iW-k+2 ■ ■ ■ Wk~\Wk and corresponds to one disjoint area in the phase 
space of a sufficiently separated disk system. If the rectangle is inside the forbidden 
region, the symbol string W-k+iW-k+2 ■ ■ -Wk-iWk is forbidden in the alphabet. By 
excluding all the forbidden symbol strings for a given number k we can redefine the 
symbolic dynamics in terms of a new alphabet with a new grammar, or in terms 

of a finite Markov graph @' 3. This is an approximation of order k to the correct 
complete alphabet describing the billiard, and this gives an approximation to a 
Markov partition of the billiard. The implementation of such approximations in 
practice will be treated elsewhere. 

3.1 Pruning for 3 disk billiard 

The pruning in the 3 disk billiard starts when the points in A with the largest value 
of |0| reach the value |0| = tt/2. These outermost points in A are the hetroclinic 
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points created by the crossing of the unstable (stable) manifold of the fixed point 
W = and the stable (unstable) manifold of the fixed point W = 1. The numerical 
value of the critical distance is r c = 2.04821419 . . .. The pruning front for 3 disks 
with r < r c is numerically determined by choosing points on disk 2 and starting 
the orbit by glancing off the disk at point x with the angle = ±7r/2. The first 
bounce in one direction is off disk 1 and in the other direction the bounce is off 
disk 3. Let Sq — 1, and a point on the first pruning front is determined by choosing 
Si = 2 and S2 = 3. The other values of s t are obtained from the numerical bouncing 
with starting point x on disk 2. To find a point on the second pruning front we 
do not include the bounce in disk 2 but assume the orbit grazes this disk. This 
gives the symbols s t = s t for t < and s t = s t +i for t > 0. The points (7, 5) are 
computed from the symbol sequences of orbits with starting point x between 27r/3 
and 7T. Plotting the points (7, 5) in the symbol plane gives fig. for parameter 
r = 2 (touching disks). The pruning front is a monotonously increasing line from 
(0.375, 0) to (0.5, 0.5) and a monotonously decreasing line from (0.5, 0.5) to (0.75, 0). 

A numerical test of this pruning front can be performed by starting one orbit at 
a random point in the billiard, letting it bounce in the billiard for some time and 
calculate the point (7, 5) for each bounce. In fig. [3] a) each bounce for an orbit with 
10 6 bounces is marked by a point in the symbol plane. As expected, the primary 
forbidden region is white in fig. |G| a) as there are no bounces with these (7, 5) values. 
In addition, there are 3 copies of this white region; the points (5, 7) that are a result 
of the time reversion symmetry (t — > —t) and the points (1— 7, 1—5) and (1—5, 1— 7) 
that are a consequence of the symmetry between odd and even in algorithm ([[]). 
This last symmetry is the results of that our symbol plane includes two fundamental 
domains. All other white areas in fig. |6| a) are images or preimages of one primary 
forbidden region, where the iteration is a shift operation on the symbol string. 

Fig. [7] and fig. |6| b) are magnifications of the symbol plane showing respectively 
the pruning front and the distribution of points of a long orbit. The long chaotic 
orbit visits points close to both pruning fronts but the probability is different as 
visualized in fig. |]b). The orbit is often close to grazing a disk but very seldom does 
the orbit bounce with angle close to n/2. The tangential bounce is very unstable 
and this gives low probability density in the symbol plane to these bounces. The 
stability is however not sensitive to how close an orbit is to graze a wall, and the 
probability density in the symbol plane approaches the second pruning front rather 
uniformly. 

One can convert the topological pruning front in fig. [7|into a symbolic description 
in different ways. The following method overcounts the number of allowed orbits. 

In fig. |7| rectangles with side length 2~ 7 are plotted together with the pruning 
front. A rectangle with side length 2~ k in the symbol plane corresponds to a symbol 
string of length 2k in the symbolic description. Each rectangle in fig. |7] is identi- 
fied with a symbol string w_ 6 w_ 5 . . . wqw^. If a rectangle is completely inside the 
primary forbidden region then the corresponding symbol string is forbidden. The 
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length 14 symbol strings . . . wqW7 that are forbidden for 3 touching disks 

are listed in column 2 in table |^. Finding the completely forbidden strings of length 
2k for all k e {1,2,.. .} gives a complete list of forbidden orbits. Table Q gives this 
list for k < 9 when r = 2. The shortest forbidden symbol string is of length 12 
{k = 6). In addition to the symbol strings w in table |2| the w' strings (obtained 
by letting — > 1 and 1 — >• 0) and the time reversed strings of w and w' have to be 
included as forbidden orbits. 

It is possible to include in the table of forbidden strings the finite symbol strings 
describing rectangles that are just partly inside the primary forbidden region. This 
gives an undercounting of the legal orbits. It is also possible to apply a combination 
of these two methods. Other strategies are using some periodic orbits close to 
the pruning front as an approximation, or to use some number of points that are 
exactly on the pruning front. It is presently not clear which approach yields the 
most convergent sequence of approximate grammars. 

The pruning front is not a continuous line in the symbol plane. The images and 
preimages of the pruning front cross the primary forbidden region and create open 
intervals in the pruning front. The result is that the pruning front itself is a Cantor 
set. The open intervals do not lead to any problems in converting the pruning front 
into a table of forbidden orbits. Any monotone curve in the interval, e.g. a straight 
line, gives a correct result. 

The pruning front for an open 3 disk system with r = 2.02 is plotted in figure [8] 
and for a billiard with overlapping disks with r = 1.97 in figure |^. As expected the 
figures show that the primary pruned region increases when the distance r decreases. 
The structure of the pruning front does not change dramatically when the distance 
changes. In the billiard with overlapping disks, r < 2, the point particle can not 
bounce infinitely many times in a corner of the billiard and this gives an upper and 
lower bound on 7 and 5: < /3 < 7, S < (1 — 0) < 1. These strips of width (3 and 
their images and preimages are also forbidden regions. Figure [10] shows a long orbit 
with r = 1.9. A large area in the lower left corner and in the upper right corner are 
pruned by these strips. 

3.2 Pruning for 4 disk billiard 

The critical parameter value where the pruning start for the 4 disk system is r c = 

2.20469453.... The pruning front is similar to the pruning front for the 3 disk 

124 

billiard. This billiard and the hyperbola billiard are discussed in detail in ref . . 

3.3 Pruning for a 6 + 1 disk system 

The billiard with 6 disks on a large circle and one disk in the center has pruned orbits 
when the distance between the disks is less than r c = 3.59148407 . . .. Algorithm (0) 
with N = 6 gives the well ordered symbols for these system. Figure [2] shows the 7 
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disks and a part of one orbit that is grazing disk 7 (the center disk) and the symbolic 
coordinate (7, 5) of this orbit is on the pruning front created by this center disk. 
The symbol plane used corresponds to the phase space of disk 1. In this system the 
disk 7 has a different phase space and symbol plane and the pruning front obtained 
is only valid for one of the disks on the large circle. There are two pruning fronts in 
the symbol plane; an orbit bouncing off disk 1 can graze either disk 7 as in figure [| 
or it can graze disk 2. The primary pruned region in the symbol plane is limited by 
these pruning fronts as shown in figure 11 for parameter r = 2.2. One orbit grazes 
both disk 2 and disk 7 and this orbit gives a point connecting the two pruning 
fronts. 

The primary pruned region is converted to a list of forbidden symbol strings 
similarly to the 3 disk billiard. Rectangles of size 5~ k correspond to symbol strings 
W-k+iW-k+2 ■ ■ ■ Wk-iWk and the rectangles in the primary pruned region correspond 
to forbidden strings. Table [3] lists all the forbidden orbits of length 2 and 4. 



4 Conclusions 

We have shown that it is possible to define a well ordered covering alphabet for 
different concave billiard systems with a finite number of concave walls. For these 
systems the primary forbidden region is a simple connected region in a symbol plane. 
The border of this region is the pruning front, and we have found this pruning front 
numerically for the 3 disk system with different parameter values. The pruning 
front is also obtained for a more complicated billiard with 7 disks. The general 
method to construct a symbolic alphabet ordered the same way as the ordering of 
orbits in the phase space is discussed. A well ordered alphabet is necessary when 
finding the pruning front because, for any other alphabet the primary pruned region 
is a very complicated region. 

We show how the topological information, the pruned region, can be converted 
into symbolic dynamics, a list of forbidden symbol strings. The tables list some of 
such forbidden strings. 

The method can easily be applied to other finite concave billiards with similar 
symmetries as those considered in this article. Infinite billiards like the Bunimovich 
stadium billiard, the Sinai billiard and the Lorentz gas systems are more complicated 
because it is difficult to obtain a well ordered covering alphabet which is not pruned. 
These systems will be treated elsewhere. A smooth potential described by the same 



symbolic dynamics as a disk billiard systemES' S may have a pruning front and 



will also be treated in a similar way. It is however problems in defining a bounce 
in the smooth potential and this is presently under investigation. 
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1. The three disks in configuration space enumerated anticlockwise. 

2. The 6 + 1 disks, enumerated anticlockwise, r = 2.2. The orbit is tangent to 
disk 7, and therefore belongs to the pruning front for the system. 

3. Show how two parallel orbits change orientation when bouncing in a concave 
wall. 

4. Show that two parallel orbits do not change orientation bouncing in a convex 
wall. 

5. The allowed, boundary and forbidden orbits: a) Two legal orbits, b) Two 
orbits on the pruning front, c) Two forbidden orbits. 

6. The symbol plane picture of an orbit bouncing 5- 10 6 times in the closed 3 disk 
system with touching disks (center distance r = 2). Each bounce corresponds 
to a point in the plane, a) the whole symbol plane, b) magnification of a 
part of the symbol plane showing the primary pruned region. 

7. A magnification of the pruning front for 3 disks with r = 2, with line mesh 
spacing at 2 -7 . The symbols give the symbolic past and the symbolic future 
of the rectangles. 

8. The pruning front for 3 disks with r = 2.02. 

9. The pruning front for 3 disks with r = 1.97 

10. An orbit bouncing 2- 10 6 times in the closed 3 disk system with center distance 
r = 1.9. 

11. The pruning front for 6+1 disks, r = 2.2. 
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S S 1 S 2 S 3 S 4 


W1W2W3W4 


12121 


0000 


12123 


0001 


12132 


0010 


12131 


0011 


12313 


0100 


12312 


0101 


12321 


0110 


12323 


0111 


13232 


1000 


13231 


1001 


13213 


1010 


13212 


1011 


13121 


1100 


13123 


1101 


13132 


1110 


13131 


1111 



Table 1: The table shows the ordering of 4 bounces of a particle starting on disc 
1 at time 0. The symbol s t is the number of the disc, while Wt is new symbols 
reflecting the ordering in phase space. 
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_U7_5 • • • U>6- 


-W-Q . . . U>7_ 


_U7_7 . . . W£ 




000000 • 100000 


0000000 


• 1000011 


00000000 ■ 


-I /~\ /~\ /~\ -I /~\ -I /~\ 

10001010 




0000000 


• 1000010 


00000000 ■ 


-I /~\ /~\ /~\ -I /~\ /~\ -I 

10001001 




0100000 


• 1000001 


00000000 • 


10001000 




0100000 


• 1000000 


01000000 • 


10000101 




1100000 


• 1000000 


01000000 • 


10000100 




0010000 


• 1000000 


11000000 • 


10000100 




1000000 


• 0111111 


01100000 ■ 


10000010 




0000000 


• 0111111 


11100000 • 


10000010 




0000000 


• 0111110 


U1U1UUUU • 


1UUUUUU1 




0100000 


• 0111111 


01010000- 


10000000 








11010000- 


10000000 








00110000- 


10000000 








10110000- 


10000000 








01110000- 


10000000 








11110000- 


10000000 








00001000 • 


10000000 








01010000 • 


01111111 








10010000 • 


01111111 








00010000 • 


01111111 








11100000- 


01111111 








01100000- 


01111111 








01000000 • 


01111101 








00000000 • 


01111011 



Table 2: The forbidden orbits in 3 discs, r = 2. "Under counting" approximation 
to level 9. 
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_W_ 8 • • • Wg. 








nnnnnnnnn 
uuuuuuuuu 


m 1 1 1 nm 1 

Ull llUUll 


1 i m ni nnn 


1UUUUUUUU 


UU1UUUUUU 


1 r\r\r\r\ 1 1 1 n 
1UUUU111U 


nnnnnnnnn 
UUUUUUUUU 


m i i i ni nn 
Ul 1 11U1UU 


ni ni ni nnn 


i nnnnnnnn 
1UUUUUUUU 


nm r\r\r\r\r\r\ 
UU1UUUUUU 


1 nnnm i 1 1 
1UUUU1111 


nnnnnnnnn 
UUUUUUUUU 


U1111U1U1 


1 nn 1 n i nnn 


i nnnnnnnn 
1UUUUUUUU 


1 i r\r\r\r\r\r\r\ 
11UUUUUUU 


1 nnm r\r\r\r\ 
1UUU1UUUU 


i nnnnnnnn 
1UUUUUUUU 


U1111U1U1 


nnn 1 ni nnn 
UUU1U1UUU 


i nnnnnnnn 
1UUUUUUUU 


m nnnnnnn 
UlUUUUUUU 


1 nnm nnnn 
1UUU1UUUU 


m nnnnnnn 
UlUUUUUUU 


m 1 1 1 m 1 1 
U1111U111 


111 nm nnn 


i nnnnnnnn 
1UUUUUUUU 


ni nnnnnnn 
UlUUUUUUU 


1 nnm nnm 
1UUU1UUU1 


nm nnnnnn 

UU1UUUUUU 


ni i i i 1 nnn 

Ull 11 1UUU 


ni i nm nnn 
U11UU1UUU 


i nnnnnnnn 
1UUUUUUUU 


UUUUUUUUU 


1 nnm ni i n 
1UUU1U11U 


nm nnnnnn 

UU1UUUUUU 


ni 1 1 1 1 nm 

Ull 11 1UU1 


i ni nm nnn 


i nnnnnnnn 
1UUUUUUUU 


UUUUUUUUU 


1 nnm m 1 1 
1UUU1U1 11 


m i nnnnnn 

Ul 1UUUUUU 


niiiiinii 

Ull 11 1U11 


nm nm nnn 
UU1UU1UUU 


i nnnnnnnn 
1UUUUUUUU 


UUUUUUUUU 


1 nnm i nnn 
1UUU11UUU 


111 nnnnnn 
111UUUUUU 


niiiiinii 

U11111U11 


1 1 nnm nnn 
11UUU1UUU 


i nnnnnnnn 
1UUUUUUUU 


UUUUUUUUU 


1 nnm i nni 
1UUU11UU1 


nnm nnnnn 

UUU1UUUUU 


niiiiinii 

Ull 11 1U11 


m nnm nnn 
U1UUU1UUU 


1 nnnnnnnn 
1UUUUUUUU 






nm i nnnnn 

UU11UUUUU 


ni 1 1 1 1 1 ni 

U111111U1 


nnn 1 1 r\r\r\r\ 
UUU11UUUU 


i nnnnnm n 
1UUUUUU1U 






i ni i nnnnn 
1U11UUUUU 


ni 1 1 1 1 1 ni 

U111111U1 


111 m r\r\r\r\ 


i nnnnnnn n 
1UUUUUU1U 






mil nnnnn 
Ul 11UUUUU 


ni 1 1 1 1 1 ni 
U111111U1 


ni i m nnnn 

U 11U1UUUU 


i nnnnnm n 

1UUUUUU1U 






1111 nnnnn 


01 1 1 1 1 1 01 

Wllllllul 


011010000 


100000011 






minmnnn 


011111110 


010010000 


100000100 






mininnnn 


011111111 


100010000 


100000100 






111010000 


011111111 


000010000 


100000100 






000110000 


011111111 


000010000 


100000101 






100110000 


011111111 


001100000 


100000110 






010110000 


011111111 


010100000 


100001000 






110110000 


011111111 


100100000 


100001000 






001110000 


011111111 


000100000 


100001000 






101110000 


011111111 


000100000 


100001001 






011110000 


011111111 


111000000 


100001010 






111110000 


011111111 


011000000 


100001010 






000001000 


011111111 


001000000 


100001100 






100001000 


011111111 


001000000 


100001101 







Table 2: Continue. 
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-WoWl- 




-i 


• • -W 2 - 










4-1 


04 


20 


43 


•20 


32 


11 


12 


14 


11 


14 


00 


14 


3-1 


14 


20 


42 


•20 


32 


12 


02 


13 


21 


14 


10 


14 




24 


20 


32 


•20 


32 


13 


02 


14 


31 


14 


20 


14 




34 


20 


42 


• 10 


32 


14 


01 


13 


41 


14 


30 


14 




44 


20 


42 


• 11 


22 


11 


11 


13 


00 


13 


40 


14 




03 


20 


42 


• 12 


22 


12 


21 


13 


10 


13 


10 


10 




13 


20 


42 


• 13 


22 


13 


31 


13 


20 


13 


00 


10 




23 


20 


42 


• 14 


22 


14 


41 


13 


30 


13 


00 


11 




33 


20 


32 


• 10 


12 


13 


01 


14 


40 


13 







Table 3: The forbidden orbits in 6+1 discs, r = 2.2. "Under counting" approxi- 
mation to level 2. 
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